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functions of the signal-to-noise and signal-to-interference 
power ratios. The results for coherent detection of MFSK 
(M=2,4,8) in the presence of additive white Gaussian noise 
show excellent agreement with the theory. On the other hand, 
simulation results for the probability of bit error for non- 
coherent detection of MFSK differ (-21% average) from the 
theory suggesting a possible ‘systematic’ error in the 
Communications Toolbox implementation of the non-coherent 


MFSK detection. 
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E. INTRODUCTION 


A. OBJECTIVES 

In many military scenarios a large number of different 
communication systems operate simultaneously from the same 
site or platform, ideally with as little mutual interference 
as possible under the constraints of available frequency 
allocations and the physical separations of the antennas. It 
is important to predict the effects of co-channel 
interference such that measures may be taken to reduce the 
interference to a tolerable level. 

The objective of the thesis is to model an M-ary 
frequency-shift keying (MFSK) communication system in the 
time domain using S jen eae tke and MATLAB/SIMULINK 


Communications Toolbox. MFSK performance is verified for 


additive white Gaussian noise (AWGN). In addition, MESK 
performance mer Simultaneous AWGN and co-channel 
interference is determined. This is accomplished by 


obtaining estimates via computer simulation of the bit error 


rates under these conditions. 


B. METHOD 
Using Simulink and the Communications Toolbox allows 


one to concentrate on system modeling rather than on the 


details of code development in Matlab or some other computer 
language. The block-diagram based SIMULINK time domain 
modeling iS a convenient tool for the visualization of 
communication Signals at various stages of transmitters and 
receivers and for “Monte Carlo" type simulations of their 
performance. The results of Monte Carlo type computer 
Simulation were used to establish the bit error rates, under 
realistic conditions of noise and interference and for 
different values of  transmitter/receiver and channel 
parameters such as the power ratios of the desired and 
interfering signals. 

In Chapter II the Simulink and Communications Toolbox 
is described. Chapter III discusses the models used for the 
simulation. The simulation results for coherent and 
noncoherent MFSK (M-2, M-4, M-8) in the presence of AWGN are 
verified in Chapter IV, and in Chapter V the results are 
discussed for the case involving the presence of both AWGN 
and interference. Finally, Chapter VI presents a summary and 


conclusions. 


II. SIMULINK AND COMMUNICATIONS TOOLBOX 


A. INTRODUCTION 

Simulink is a program for modeling and simulating 
dynamic systems including linear, nonlinear, discrete-time, 
continuous-time, and hybrid systems. As an extension to 
Matlab, this environment adds many features specific to 
dynamic systems while retaining all of Matlab’s general- 
purpose functionality. Model definition and model analysis 
are the two steps in a Simulink computer simulation. In 
practice, these two steps are often performed interactively 
as the model is created and modified in order to achieve the 
desired behavior. 

Simulink represents systems as block diagrams. Using 
mouse driven commands, one can create models whose 
parameters can be edited or changed using the keyboard. 
Also, operations can be clicked-and-dragged, and results can 
be displayed in real time during a simulation. The models 
can be analyzed either by choosing options from the Simulink 
menus or by entering commands in Matlab’s command window. 
The results, which can be viewed while the simulation is 
running, are made available in the Matlab workspace after a 
Simulation has been completed. Since Simulink is built on 


top of the Matlab environment, it has many capabilities, 


which can be further enhanced with application-specific 
POOrsoxXKes = Generally Speaking, Simulink combines the power 
and ease-of-use of an application package with the 


flexibility and extensibility of a language [Ref 2,1]. 


B. TIME DOMAIN SIMULATION 

In order to evaluate the time domain simulation, the 
terms “modules” and “single-point format” need to be 
explained first. The modules of a simulation program can be 
regarded as subroutines that operate on the data, which is 
provided as input. The program calls the first module, which 
processes the data and returns the results. Then it calls 
the second module to process the data delivered from the 
first module and so on. This continues until all modules 
have been called, upon which the program repeats from the 
beginning with new input data. The term “single-point” 
refers to the amount of data that can be processed on each 
call of a module. When only one data point (only one time 
instant) can be processed or one data point is passed for 
computation from module to module, such simulation format is 
called “single-point data.” On the other hand, if a whole 
block of N data points (N time instants) is processed on 
each call, the format is called “block data.” This means 
that for an equal amount of total data, single-point format 


calculations are substantially more time consuming than the 


block data operations since single point data calculations 
incur at least N times as much calling overhead as block 
data simulations for the same (N) number of time instants. 
The single-point format (which is the only alternative when 
a System includes feedback loops) processes only time-domain 
Signals (periodic or aperiodic). Usually, single-point data 
is regarded as one sample from a long, aperiodic signal, and 
thus a single-point simulation will exhibit the same 
transient behavior as the system being simulated. Also, 
Single-point simulations can be arbitrarily long in duration 
without special provisions. This makes the simulation more 
convenient when collecting data for performance statistics 


[Ref 6]. 


e CONSTRUCTING A SIMPLE MODEL 

Simulink uses block diagrams to represent dynamic 
systems. Defining a system is much like drawing a block 
diagram. The “building blocks” can be copied from various 
block libraries, which can be either built-in (Simulink or 
Communications Toolbox libraries) or user-created. The 
built-in Simulink is organized into sub-libraries according 
to block function or nature. In order to assemble a system, 
blocks can be copied from an existing library and connected 
using the mouse and new block parameter values inserted in 


place of the default values, thus saving time when building 


new models. As an example, the process of setting up a 
simulation of a very simple system is outlined. Suppose that 
the output of a pulse generator is to be observed on an 
oscilloscope. First, the command Simulink is typed at the 
Matlab command window prompt in order to start Simulink and 
access the built-in Simulink libraries. A double-click on 
these basic libraries opens another window containing the 
related blocks. Figure 1 shows the basic Simulink libraries 
(whrchweshow upon. typing Sethe command Simudank). «Figure 2 


shows some of the blocks available in the Sources library. 
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Figure 1. The Standard Simulink Block Library. 
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Figure 2. Blocks Available in the Sources Block Library. 


TO build the simple simulation the Sources sub-library 
needs to be opened. This displays yet another Simulink 
window from which the discrete pulse generator can be copied 
to a “New File” window. Similarly, by opening the Sinks sub- 
library, the oscilloscope block can be copied and placed 
into the new system window. Figure 3 shows the discrete 


pulse generator placed into a new window. Figure 4 shows the 


oscilloscope block placed into the new window and finally 
Figure 5 shows the finished model, with the source connected 


to the oscilloscope. 
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Figure 3. Placing the Discrete Pulse Generator Block. 
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Figure 4. Placing the Oscilloscope Block. 
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Figure 5. Completed Simulink Model. 


All copied blocks maintain the original default 
internal parameters until modified by the user. Therefore, 
by simply connecting the sine wave "source" block to the 
oscilloscope "sink" block and selecting the parameters of 
the source, the system is ready to perform a simulation [Ref 


ma 


D. COMMUNICATIONS TOOLBOX 

The Communications Toolbox is designed for use with 
Matlab and Simulink and represents a collection of 
computational functions and simulation blocks for research, 
development, system design analysis, and simulation in the 
communications area. The toolbox contains ready-to-use 
functions and blocks, which can be easily modified by the 
user to implement the required schemes, methods, and 
algorithms. Simulink blocks and Matlab functions accelerate 
the design process by helping users to rapidly develop and 
analyze different system designs. By using the 
Communications Toolbox (whose functions can be either called 
directly from the Matlab workspace and .m files or 
implemented as specialized Simulink blocks) time-domain 
Simulation of various communication systems can be 
implemented. The toolbox adds a variety of Matlab functions 
and Simulink blocks that are very useful for communication 


system simulation. These include: 
e Signal Generators (data source) 
e Source Coding/Decoding 
Eco Control Coding 
e Baseband and Passband Modulation/Demodulation 


e Transmitting and Receiving Filters 
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e Baseband and Passband Channel Models 


Multiple Access 
e Syncaronizacion 


e Galois fields calculations. [Ref 1] 


E. CONSTRUCTING A COMMUNICATIONS SYSTEM MODEL 
A simple communication system contains these basic 

components: 

e information source, 

e the transmitter, 

e the channel, 

e the receiver, and 

6 information “sink.” 
Suppose that the performance Of a coherent MFSK 
communication in the presence of additive white Gaussian 
noise needs to be determined. The Communications Toolbox is 
started by typing the command commlib at the Matlab command 
line prompt. This gives access to the "top level” 
communications library shown in Figure 6. The transmiting 
part of a communication system is comprised by the upper row 
of blocks while the receiving part is comprised of the 


FOottom row ©f bloeks- 
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Figure 6. The Communications Toolbox Simulink Block Library. 


Since in most cases the receiving computation is 
exactly the reverse of the transmitting computation, a 
double click on any of the blocks in this library opens a 
sub-library that includes both the transmitting and the 
corresponding receiving blocks. A sub-library typically 
contains function blocks used to build the models and demo 
blocks (cyan colored) that are used to demonstrate a block 
operation. Finally, by opening the source/sink library 


various types of sources become available. In order to build 


I 


A 
aif Bone m, d f » "m cm 


an MFSK system the "Sampled Read from Workspace" block is 
copied from the Source sub-libarary and placed into a new 
Simulink window. The "Sampled Read from Workspace" block 
reads a row of data from the Matlab workspace at every 
sampling point. Similarly, by opening the  modulation/demo- 
dulation library (Figure 7) and double clicking mthe 
Digital Modulation/ Demodulation sub-library, the MFSK Mod 
Baseband and the Coherent MFSK Demod blocks can also be 


copied and placed into the new window. 
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Figure 7. The Modulation/Demodulation Library. 
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The Coherent MFSK Demod block demodulates the input, 
which is a complex analog signal (because of the baseband 
implementation of the passband modulation) and is defined by 
five parameters:[M-ary Number, Tone Space, Symbol Interval, 
Initial Phase, Sample Time]. The MFSK Baseband Mod block 
accepts a scalar input in the range [O0 M-1] and is defined 
by four parameters: [Tone space, symbol interval, initial 
phase, sample time]. 

From the channel sub-library the AWGN channel block is 
selected and copied into the new window. This block adds 
AWGN to the signal and is defined by three parameters. The 
first two specify the mean and the variance of the noise 
output and the third initializes the seed for the Gaussian 
random number generator. Finally, the error-rate meter block 
is copied from the sink sub-library. This block, shown in 
Figure 8, performs the symbol-to-symbol comparison between 
the sender and the receiver, counts symbol errors, and 
computes the symbol error rate. As shown in Figure 8, 50 
symbols have been transmitted and the receiver has detected 


13 symbols incorrectly. So the symbol  error-rate is 
l 
= = 0.270833 . All the copied blocks are connected in order to 


complete the MFSK system shown in Figure 9. 
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Figure 8. The Error-Rate-Meter block. 
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Figure 9. A Simple Model of a Coherent MFSK System. 
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E. MODEL ANALYSIS 

Comunication system simulation can take place upon the 
completion of the communication system model. There are 
basically three different ways of implementing and 
controllindgENassumiilucss cw hc fdrstedsetosüm the simulation 
interactively from the menu bar and observe the signals at 
various inputs and outputs using oscilloscopes while the 
simulation is running. This method is easy to learn and 
simple to use, but it also incurs the highest overhead and 
is therefore the slowest to execute. The second option is to 
run the simulation by executing a .m file. The main 
advantage of this method is that it provides greater 
flexibility than the first method. For example, parameters 
in the blocks can be changed for repeated simulations and 
initial conditions for various blocks can be overridden. 
Running a simulation from a .m file allows one to change one 
parameter (or parameters) each time the simulation is 
repeated and the results of the simulation can be brought 
into Matlab workspace for further analysis. Finally, the 
third possibility of implementing a simulation (which is the 
most complex and flexibile) is based on the fact that every 
Simulink block has a Matlab equivalent as an S-function. 
Each such S-function has the same name as the corresponding 


model and can be called in a variety of ways to provide 
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information about the model it represents. The S-functions 
can provide information about the number of inputs, outputs, 
and states (both continuous and discrete) of the model as 
well as the derivatives of the outputs. All of the analysis 
tools provided with Simúlink interact witi modelbmihr ough S- 
functions. The linmod and trim .m files are examples, which 


use these functions. These are shown in Table 1. 


i linmod Generate linearized model of a 
continuous system 
2 


SERT Generate linearized model of a 
B a - with discrete elements. 
3 tr im Trim systems about equilibrium 
ae 


Table 1. Partial List of Analysis Functions (from Ref [9]). 





Generally speaking, the three ways to implement/run a 
simulation are not strictly delineated and different usage 
may be appropriate at different stages of model development 
(Ret eZ, 6, Iie 

Finally, Simulink provides a number of methods for the 


numerical solution of ordinary differential equations that 


1 


can be selected prior to simulation and changed even while 
the simulation is running. These numerical methods for 


solving "initial value problems" are shown in Table 2. 


ls am Linear systems method 
2 Y 2" order Runge-Kutta method 
AS 4^ order Runge-Kutta method 


euler Euler's method 


Adams predictor-corrector method 
gear Gear's predictor-corrector method 
for stiff systems 


Table 2 maune ei ons used for Simulation of Ordinary Dif 





ferential Equations (from Ref. [9]). 
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III. MODEL OF A MFSK PASSBAND DIGITAL COMMUNICATION SYSTEM 
AND ITS BASEBAND EQUIVALENT 


A. BASEBAND EQUIVALENCE FOR PASSBAND SIGNALS AND SYSTEMS 

The most convenient way of representing a narrowband 
system for simulation purposes is to work with equivalent 
baseband quantities. Consider a communication system 
operating with 1 GHz andja bandwidth of 20 MHz. According to 
the Nyquist Sampling Theorem, we need a sampling rate of at 
least 2Xfmax = 2x1l,010 MHz to simulate this system. A 
sampling rate of this magnitude is not practical. Since the 
bandwidth of the baseband equivalent is several orders of 
magnitude smaller than the passband  signal's center 
frequency (20 MHz versus 1 GHz), we can simulate the 
passband system by establishing its baseband equivalent 
system and simulating the equivalent system using a 20 
Megasamples/second sampling rate instead of the over 2 
Gigasamples/second sampling rate required for the passband 
Signal. Therefore, simulating the baseband equivalent of a 
passband signal substantially reduces the required computer 
storage (memory and disk space) and the calculation time. 

A (real) passband signal can be expressed in polar 


(amplitude and phase) form as: 
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pc sosI2mTf.t ^ Q(t)1 => 


Reta (O ed" ¿trote Y M 


> 
Chi 
li 


unes alccnhe.envelope,j f. isothe carrier frequency, Re 
is the real operator and g(t) is the phase of x(t). 


A (real) passband signal can be also expressed in 


Cartesian form: 


Pee) cr (EOS (2 UE kalt) ssn(2m-tf.t (Seal) 


where X(t) is the in-phase component and x(t) is the 
Quandrature component of x(t). 
The corresponding complex baseband (or lowpass) signal 


is defined as: 


p M ceci : (35938 
CUTS 


eee) T E E) (3.2b) 


In order to transform the real passband signal of 


equation (3.1a) to the complex baseband signal of equation 


(Gelb yeversuppress thevcarcricr by Gropping ©. c. 
The term a(t)e?'') is called the complex envelope of 


x(t). The baseband nature of X(t) and the low sampling 


rate required to satisfy Nyquist's Theorem for narrowband 
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signals is evident. Complex envelope analysis is very useful 
in analyzing the responses of passband systems to passband 
signals. Finding the response of passband systems to 
passband signals is one of the most essential tasks in 
communications system analysis. 


If we denote the Fourier Transform of the passband 
Signal x(t) as X(@) w=2nf and the Fourier Transform of its 
complex envelope x.(t) as X.(0) then: 

X(0)=0.5 X.(@-@) + 0.5 X.* (-@-@,) (Soa) 
where * represents the complex congugate. The first. term 
represents complex envelope’s baseband spectrum scaled by 


0.5 and shifted to the carrier frequency. The second term 


represents the complex envelope's baseband spectrum scaled 
by 0.5, shifted to carrier frequency, folded about o-0 axis 


(to negative frequencies) and conjugated. Also, we can 


conclude that the first term represents the spectrum of 
X(@) for positive frequencies while the second term 


represents the spectrum of X(@) for negative frequencies. 


Baseband equivalents of passband systems can be 


determined in a similar manner. For example, consider a 
filter with a passband frequency response H(@) and real 


impulse response h(t). We can now express the frequency 


Za 


response H(0) in terms of the baseband equivalent frequency 


response: 


H(@)=0.5 H.(@-@,) + 0.5 H.* (-@-@,) (3.3D) 


hemi egua riena e a a. 3b); we Can Conclude that 
the complex envelope of the passband response to a passband 
input can be obtained as the response of an equivalent 
baseband system to the complex envelope of the passband 
input x(t). In this manner, the baseband equivalent system 
can be used to avoid calculations involving high frequency 


carriers [Ref 5,8]. 


B: MONTE CARLO SIMULATION 

Symbol error probabilities (SEP) and bit error 
probabilities (BER) can be estimated using simulations with 
very large numbers of symbols transmitted by counting the 
symbol errors and taking the ratio of the number of symbol 
errors to the number of symbols as the number of symbols 
tends to infinity, assuming that the symbols are generated 
independently of one another. This is the essence of the 
Monte Carlo Simulation technique. The technique is simple 
but requires significant computing power and/or long 


simulation times when accurate estimates are needed for 
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small symbol/bit error rates. For example, if we count 10 
errors, we can find -from Figure 10 the SER intomo ls HE 
contain the true SEP with 90$, 95$, and 99$ confidence {Ref 


om 


po» 


E ue. 


= 10” 90% 95% 99% 


ZZ eo 
] get) 
10" po 10 108 
Number of symbols 


Figure 10. Confidence bands on P, when observed value is 10" 


for the Monte Carlo technique (source: Jeruchim etal.,1984). 


Erom Figure 10, 1f a symbol (or Dit) error probability 
MSO yl) en | form aloe experimental errors the actual 


(“true”) probability of error is between 0.1 to 0.065 and 
E O O aV Oaa ean l den e [However Lr errors are not 


independent from each other (errors occur in a "burst"), the 
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Sierro: Im a burst is more informative of the error 
statistics than all subsequent errors in the same burst. 
Therefore, in the “dependent error” case we need to count a 
larger number of errors for the same SEP confidence levels, 


which means even longer simulation times [Ref. 5]. If the 
number of experimental errors is denoted Me and the number 


of the received bits (or symbols) is denoted N,, then the 


confidence with which the BER is upper and lower bounded by 


Nu "n ; 
N. and NO! respectively, is given by 
b b 


n Nu 
Prob [—- «x BER x —/ — (9 
ES m n.j (314) 


b b 
where ©(0<0<1) is the confidence level and ft, Mu is the 
confidence interval. For a given number of experimental 
errors and if certain reasonable conditions are fulfilled 
(effects of disturbances are independent for each error and 
error-count statistics can be approximated as Gaussian 
distributed), then the lower and upper number bounds for the 
errors with 95.4% confidence interval are given from the 


following equations: 


Nu=Ne +2-2n, +1 (SON) 
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Nu=Ne +2+2 yn, +1 (3x26) 


where Ne 1S the number of experimental errors and Nu is 


the confidence interval. Generally speaking, error bounds 
converge slowly to each other (that is the estimates become 


more accurate) as the number of observed experimental errors 
Ne increases. For ten experimental errors, the true number 


of errors is between 5.4 and 18.6 (between approximately - 
50$ and +76% of the observed value), with 95.4% confidence, 
while for 100 experimental errors the true results will be 
between 81.9 and 122 (between approximately -18% and +22% of 


the observed value) also with 95.4% confidence. The values 


of n and Mu from (3.5) and (3.6), respectively, and their 
ratio are tabulated in Table 3 for several values of fh and 


mer confidence” level © = 907954 (95.45 contidence) [Ref sol. 


zu 


Experimental | Lower limit Upper limit of 


errors Cf Venue ter= ~rue errors (Ny) 


rors (N) 





Table 3. 95.4% Confidence Intervals for Various Number of 
Experimental Errors (from Ref. [6]). 


C. M-ARY FREQUENCY SHIFT-KEYING (MFSK) 

M-ary frequency-shift keying (MFSK) modulation does not 
require phase coherent detection. This means that MFSK can 
be used in applications where the phase shift incurred by 
the transmission channel varies rapidly, as is the case, for 
example, in communications with high-speed aircraft. With 
MFSK modulation, each symbol is assigned a specific 
frequency and an MFSK signal can be described by the 


following equation: 


p mE EE EnpPeessT2t toan) m cR 
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where: 


Af is the spacing between any two adjacent frequencies: 
Af=f.-Éan-1 (mia) 


ay; is determined by the information sequence and takes M 
integer values 
fo is the lowest frequency 
A is the amplitude of the MFSK signal 
Pa is the phase angle 
p(t-nT) is the unit-amplitude rectangular pulse. 

qe Bit Error Probability For Coherent MFSK 

In order to obtain the exact expression for the MFSK 
bit error probability, the exact expression for the MFSK 
symbol error probability must be first derived. The coherent 
MFSK symbol error probability is given by the following 


expression [Ref 7]: 


» [2SNR sym 


sym = ER 5 seu): 
PREIS = laca Wo 72) n dc (3.8) 


where: 


2 


C-z/o (oa 
O^ is the variance of wideband Gaussian noise 


M is the number of symbols 


SNR. 1S the symbol SNR 





log(M) 
SNR — SNRyp;.———— 3.8D 
sym DIt log(2) ( ) 
ando) is the O-fúnction given by the following: 
Q (x) = La EI = m (ION 
2 qoi 
where the error function (erf) is defined as 
1 ^ m 
eriz) = e dx (38a 


Jas 


The following auxiliary functions are introduced for 


convenience: 
| log( M) 
= 2 — A SNR. : 
MS SNE»b:c) log(2) bit (3 9) 
i ~(6-y(M,SNR,,))? 
w(C,M,SNRpit) = "- WA —— (3. 10) 
JE 


Z8 


By using these auxiliary functions we now can express the 
bit error probability [Ref 7] as a function of the number of 


symbols M and the bit signal-to-noise ratio SNRpi:: 


Porc 2125 M N M,SNR 3.11 
s ses VIS] Y(M,SNR,, )-5 WG vit AG C 
2 Bit Error Probability For Noncoherent MFSK 


In order to obtain the exact expression for the MFSK 
bit error probability, the exact expression for the MFSK 
symbol error probability must be first derived. The symbol 
error probability for noncoherent MFSK is given by the 


following expression [Ref7]: 


=A 


€ CD*M-D' uix, 


ic (3.12) 
£4 (k 4- D)k'(M -1- k)! 


Pe (sym) = 


where: 


M is the number of symbols, M = E aa ene 


N , ! 
= is a two-sided power spectral density,  No- 0/4T, and 


A is the amplitude of the MFSK signal. 
We can also express the noncoherent MFSK symbol error 


probability in terms of the symbol energy: 


DI 


ELE 
EDT MN C 
= CD*(M-D! m m. 


M =} (Kk - DK(M -1- Kk)! 
or in terms of the symbol signal-to-noise ratio: 
M1 (-DF(M-D!  -—sm,, 
Pe (sym) = ee a (SAN 


E (k+DKIM-1=k)! 


The symbol error probability for noncoherent MFSK is 
related to the bit-error probability for M = 2" by [Ref7]: 


-k log M) 








M-1 — K — | "me! JSNR bn: 
Po (MM, SN) = TEMS js DM =D! a log(2) ese 158 
2M-1 o (K +1)k!(M -1-k)! 
OF 
-D*(M-D! a 
Pewoit) ( M,;SNRpit) = iu zx Y TD e Dc ] (3.59) 
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De MFSK MODELS AND BLOCK ANALYSIS 
a Coherent MFSK Model (Noise Only) 
The following figure shows the coherent MFSK (noise 


case only) model, which has been used for the simulation. 
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Figure 11. Coherent MFSK Model for Noise Only Case. 


2. Block Analysis 


a) Coherent MFSK Mod Baseband 

Four parameters define the coherent MFSK Mod 
Baseband block: 
[Tone space, symbol interval, initial phase, sample time]. 
The MFSK Baseband Mod block accepts a scalar input in the 
range [0, M-i]. (For example, M in the coherent BFSK model 
357). <Ihe™otrput of the MFSK Mod block is a Unit-amplPeme 
FSK modulated analog signal. The following must be specified 


for the MFSK Mod block: 
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e Tone space (Ag) (also known as the frequency 
separation between two neighboring frequencies and 
defined as  A;1/(2Ts), where T, is the symbol 
interval). 

e Symbol interval T, (also the data sampling time) 

e Initial phase (the initial phase of the carrier, 
commonly set to 0) 

e Sample frequency (the inverse of the signal sampling 
interval) 


The theoretical BER for coherent BFSK is given as: 


n-d e| (3.17) 
O 


where Eg is the average bit energy and No is the single- 
sided noise power spectral density. 

From the Nyquist Theorem, at least two samples per 
shortest period are needed in order to represent a signal. 
For the baseband equivalent model of binary FSK (BFSK), one 


symbol is represented as DC and the other symbol as a 
carrier with the frequency Af. We have selected the symbol 


interval T, as 1 second. (In the binary case this is also 
the bit interval Ty). For a good representation of the 


Signal we sample at twice the Nyquist rate; that is, four 


eZ 


samples per bit interval, or four samples per second, are 


provided. The number of samples per second N, is given as: 


T, 
N= — PEPE, 
A, 
where A, is the MFSK signal (not data) sampling time. For 
BFSK and sampling at twice the Nyquist rate, we get: 
AG 
A.= + ERIE 


Or 


^= lsec/4samples per sec > A, = 0.25Hz. 


For orthogonal MFSK and coherent detection, the frequency 
separation must be an integer multiple of T,/2. The lowest 


separation for orthogonality with coherent detection is: 


As = x => Ac -> Tk => dems => Ag = 0.5Hz (320) 


This value is defined as the “tone space” (in Hz) in the 


Simulink model. In the case of a coherent 4FSK or 8FSK 


signal the sample time At must be changed because the 


highest frequency in the spectrum is approximately (M-1)Af 
Since the maximum frequency for MFSK is defined as fax = (M- 


1) Ag. In order to have a good representation of the signal, 
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we transmit four samples per symbol interval (or four 
samples per second since we have selected T;,-1 second), and 


the sampling frequency for MFSK is: 
Í sampling ES 4 (M-1)Af R2 


The sampling time At is the reciprocal of the sampling 
frequency: 


Anc M Me EE (3222) 


b) Coherent MFSK Dem Baseband 

This block demodulates the input, which is a 
modulated complex analog signal. Its output is an integer in 
the range [0,M-1]. Five parameters define the Coherent MFSK 
Dem Baseband block: 
e [M-ary number, tone space, symbol interval, initial 

phase, sample time] 

The above parameters should match the ones used in the 
corresponding Coherent MFSK Mod Baseband block. We have 


selected the following values: 
l 
e Tone space (Ag): ar Ts is the symbol interval 


e Symbol interval (T;,): 1 sec. 


e Initial phase: 0 
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e Sample time : 1/(4(M-1)Af) 


e M-ary number: 2, 4, or 8 


c) AWGN Channel 

This block adds AWGN to the signal being 
transmitted through this channel. Three parameters define 
the AWGN channel. The first two specify the mean and the 
variance of the noise output. The third one initializes the 
seed. The seed effectively selects a pseudorandom eee en oe 
to be generated (for the same value of the seed the same 


sequence is always generated). 


d) Symbol Counter 

The counter block counts the number of symbols. 
There are three input ports to the counter. The first port 
inputs a constant value of 1. The second port inputs the 
clock signal, which has one pulse per each symbol generated, 
and the third input port is unused ("grounded"). The value 
of 0.5 is set as threshold. Each time the clock pulse 
crosses the threshold going from O to 1 (on the rising edge 
of the clock pulse) the counter increments its output by 1l, 


thus effectively counting the input symbols. 
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e) Error Counter 

Remoto compares the delayed "original^ signal 
(first input) and the demodulated signal (second input) at 
the end of each symbol interval. The timing is controlled 
by the third input when the rising edge of the square wave 
(from the pulse generator) crosses the threshold. The fourth 
JM MEN ne ue crounded ). The output of this block 
increases by one each time the absolute value of the 
difference between the inputs 1 and 2 exceeds the specified 
tolerance, that is each time the detected symbol differs 
from the original (sent) symbol. Therefore, by setting the 
tolerance of the input difference to 0 (or a very small non- 


zero value), the symbol errors are counted. 


£) Triggered Write To File 

Four parameters define the Triggered Write to File 
block: 
[File name, data type, number of trigger pulses between 
saved data, threshold in detecting trigger signal] 
This block writes a record to a file only at the rising edge 
of the signal coming from the pulse generator (the block can 
also isolate and write all the errors to a file in case of 
burst errors). The file can be an ASCII text file or a data 
type ("integer," "float"). The third parameter (number of 


trigger pulses between saved data) determines the number of 
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rising edges that should be received between file writes. 
Setting this field to zero forces the block to write a 
record at every rising edge of the trigger. This block has 
two input ports. The first input is the error counter value 
that is recorded into a file and the second input input is 


the "write trigger" signal from the pulse generator. 


g) Triggered Write To Workspace: 

Six parameters define the Triggered Write to 
Workspace block: 
[Workspace variable, data type, number of trigger pulses 
between saved data, maximum row number, section to keep if 
overflow occurs, threshold in detecting trigger signal] 
The saved workspace variable can be either a column vector 
when the first input port (message signal) has a scalar 
input or a matrix when the first input has a vector input. 
The first element of the input signal vector is saved in the 
first column; the second element of the input signal vector 
is saved in the second column, and so on. The block can save 
the input as string variables or data. The maximum row of 
the output variable is limited by a pre-defined number, 
which can be changed during simulation. After the limit is 
reached, the block keeps the portion of the data defined by 
the "overflow" entry. Data may be written into the workspace 


at every trigger or a number of triggers may be skipped 
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depending on the value of the 


between sucessive writes, 


corresponding rpentry (zero forces a write at every trigger) 


h) K-Step Delay 


This block implements a delay that is equal to an 


integer multiple of the sampling interval that is in effect 


for the particular simulation. 


i) Sum 


This block adds or subtracts the two inputs 


depending ONR Ne Sigugselection forgthe Inputs: tio m A 


Since the objective is to compare the original and the 


detected symbols, either +- or -+ can be used. 


j) Abs 
This block's output is the absolute value of the 


This block creates a positive value of 1 each time 


suut. 
regardless of 


the original and the detected symbols differ, 


the sign (polarity) of the difference. 


k) Pulse Generator 


Four parameters define the pulse generator block: 
of period), 


[Pulse period (symbol period), duty cycle (% 


amplitude, start time] 
This block generates a unit-amplitude square wave with the 


period equal to the symbol interval. The pulse generator 
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Starts at the end of the first symbol interval such that the 
rising edges of the square wave coincide with the end of 
symbol times. Therefore, since there is a rising edge of the 
square wave at the end of each symbol interval, there is 


always a rising edge coincident with a symbol decision time. 


l) Switch 

The switch block has three input ports. The first 
port inputs the clock signal and the second input port 
inputs the absolute value of the sum of the delayed and 


demodulated signal. The third input is unused ("grounded"). 


m) Sampled Read From Workspace 

The Sampled Read From Workspace block reads a 
matrix from MATLAB workspace. In our case the matrix is a 
matrix of pseudo-random integers generated using MATLAB’s 
random integer function randint. There are four argument 
fields, delineated by commas, in the specification of the 
randint function. The first two define the size (rows and 
columns) of the matrix to be generated. The third field 
accepts either a single integer (the number of symbols M) or 
two-integer elements array [0, M-1]. In either case the 
range of random integers (the elements of the random matrix) 


is from 0 to Mel and all M values are equally probable = The 
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last field is the seed for the pseudorandom matrix (for the 
same value of the seed the same matrix is always repeated). 
Randint [number of rows, number of columns, [0 number of 
symbols -1], seed] 

ue 

Randint [number of rows, number of columns, number of 
symbols, seed] 

In the case of noncoherent MFSK (with noise only) the 
same model as for the coherent MFSK (noise case only) is 
used with only one difference. The coherent MFSK Mod 
Baseband and coherent MFSK Dem Baseband blocks are replaced 
with the noncoherent MFSK Mod Baseband and noncoherent MFSK 
Dem Baseband blocks. These blocks are defined by the same 
parameters as in the coherent MFSK Mod Baseband and coherent 


MFSK Dem Baseband blocks. Also, the frequency separation is 


Nee@or | /toawinstead of 1/ (27.5). 


K Coherent MFSK Model For MFSK With Interference And 
Additive Noise 


The coherent MFSK model with interference AWGN is shown 


OE Qu Ie 
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Figure 12. Model for Coherent MFSK with Interference and 


Additive Gaussian Noise. 


The difference between this model and the model for 
coherent MFSK with noise only is that the white blocks 
representing the interference signal have been added. The 
interference signal is assumed synchronized with the desired 
signal. That is, symbol transitions of the interference 
occur at the same instants in time as the symbol transitions 
of the desired signal. For noncoherent | MFSK with 
interference and AWGN, the same model is used (Figure 13) 


with the coherent MFSK Mod Baseband and coherent MFSK Dem 


41 










MFSK mod 
baseband channel 


Sampled read 
from wisp 








Noncoh MFSK 
demod baseband 






Triggered waite 
to misp 


D E A dd 
SEA j IE 
A ER TA 


K-step 


delay 


Pulse 
Generator 


—LIOMAAN 


Sampled read MFSK mod 
from nispi baseband1 







Tnggered write 
to nispi 





Gain 


Figure 13. Model for Noncoherent MFSK with Interference 


and Additive Gaussian Noise. 


Baseband blocks replaced with the noncoherent MFSK Mod 
Baseband and noncoherent MFSK Dem Baseband blocks and the 
frequency separation increased to 1/T.. In order to run 
multiple simulations for a number of different values of 
parameters such as the signal-to-noise ratio (SNR), Matlab 
.m files were used. For coherent MFSK, the Matlab programs 
sn preplem and sun fski.m (listed in the Appendix) were 
used. These programs can also be combined into a single 


program. For noncoherent MFSK, essentially the same Matlab 
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.m files were used, with the appropriate variables changed. 
The .m files are used to complete the three basic tasks: 
data generation and input, simulation, and recording the 


output into files, as shown in Figure 14. 


DATA 
GENERATION 
AND INPOT 


SIMULINK 
MODEL 





ui 
FIEES 


Figure 14. Flow Chart for Matlab .m File Simulation 


Implementation. 


The first task involves the choices concerning the case 
we want to test (noise or noise and interference), the type 
of MFSK (2FSK,4FSK, 8FSK...), and the simulation parameteres 
such as symbol duration T,, oversampling factor, minimum and 
maximum signal-to-noise and signal-to-interference ratios, 
mimimum number of errors acceptable, the factor multiplying 


the number of errors (to select the number of symbols to 
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try), the maximum size of the random integer array, and the 
name of the file for saved data. The second task is the 
actual simulation started by the appropriate Matlab command 
and using the selected Simulink model (coherent MFSK or 
noncoherent MFSK). The third task is recording of the 
results: the signal-to-noise and  signal-to-interference 
ratios, the number of errors, the number of symbols, and 


symbol error probability, all as ASCII data files. 
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IV. SIMULATION, ANALYSIS, AND PERFORMANCE VERIFICATION FOR 
AWGN 

A. PERFORMANCE VERIFICATION FOR AWGN 

In this chapter computer simulation results are 
presented in order to verify the performance of  MESK 
communication systems in the presence of AWGN. These results 
are the average of three computer simulation runs for each 
signal-to-noise (SNR) ratio. The value of two has been used 
as an oversampling factor (sampling rate of twice the 
minimum value). The simulations ran until at least 100 
errors were observed. The bit error probability intervals 
for the 90$, 95$, and 99% confidence levels may be found 
trom wEyguree 0 (Chapecer Dif)! = =for a ā given number of 
experimental errors and for the confidence level of 95.4%, 
the lower and upper bounds on the number of errors can be 
found from Table 3-1. The observed error number is a random 
variable, subjected to statistical scatter arising from 
different realizations (in general, a different number of 
errors is observed each time a simulation is run with a 
different set of seeds for data and noise). In order to 
prevent ‘out of memory’ errors, the data sequences were 


limited to 10° symbols for each simulation. If less than 100 
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errors were observed in a simulation, then the sequence was 
repeated until a sufficient number of errors occurred. 

fu Results For Coherent Binary FSK (BFSK) 

The theoretical and the experimental bit error 
probability for coherently detected binary FSK (BFSK) as a 


function of the SNR in dB is shown in Figure 15. 
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Figure 15. Probability of Bit Error (Theory and Simulation) 


for Coherent BFSK Detection. 
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The relative difference between the theoretical 
probability of error for coherent BFSK ande upBSb HANE 
of error estimates obtained by the simulation are shown in 


Figurer he: 


Percent Error 
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Figure 16. Percent difference between Theoretical and 


Simulation Results for Coherent Detection of BFSK. 


The convergence of the estimate of the bit error 
probability as a function of the number of symbols 


transmitted in the course of the simulation is shown in 
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Figure 17. This allows us to estimate the minimum number of 
symbols required for a particular simulation and to verify 
that minimum has been reached. For etas particular 


simulation, convergence is achieved at approximately 1.210? 


symbols. 


Convergence 


1.86654210 = 


error bit 
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Figure 17. Convergence Plot for Coherent Detection of BFSK. 
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Ihe accuracy of the simulation estimate for the bit error 
probability for coherent BFSK is shown in Figure 18. The 
convergence to zero as the number of symbols increases is 


evident. 


Convergence Accuracy Verification 
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Figure 18. Percent Difference Between Theory and Simulation 


for Coherent BFSK as a function of the Number of Symbols. 


The differences between the theoretical bit error 
probabilities for coherent 2FSK and the simulation estimates 
for various SNR’s are plotted in Figure 19 as a histogram. 


All estimates are within the -20% to +10% interval about the 
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Euoncircolvalues, with most differences very. close to 0 


(approximately +2% off). 
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Bpuxogure l9. Histogram Of Simulation Brreonms for Coherent Beem 


2. Results For Coherent 4FSK 
The theoretical and the experimental bit error 
probability for coherently detected 4FSK as a function of 


the SNR in dB is shown in Figure 20. 
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Figure 20. Probability of Bit Error (Theory and Simulation) 


for Coherent 4FSK Detection. 


The relative difference between the theoretical 
probability of error for coherent 4FSK and the probability 
of error estimates obtained by the simulation is shown in 


Figure 21. 
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Figure 21. Percent difference between Theoretical and 


Simulation Results for Moherent Detection of 4FSk. 


The convergence of the estimate of the bit error 
probability as a function of the number of symbols 
transmitted in the course of the simulation is shown in 
Figure 22. This allows us to estimate the minimum number of 
symbols required for a particular simulation and to verify 
that the minimum has been reached. For this particular 
simulation, the convergence is evident at approximately 


4.5°10* symbols. 
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Braure 22. Convergence Plot for Coherent ARSK. 


The accuracy of the simulation estimate for the bit error 


probability for coherent 4FSK is shown in Figure 23. The 


‘oscillatory’ convergence to zero as the number of symbols 


increases 1s evident. 
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Figure 23. Percent Difference Between Theory and 
Simulation for Coherent 4FSK as a function of the 


Number of Symbols. 


me differences between the "theoretical Dit erron 


probabilities for coherent 4FSK and the simulation estimates 


for various SNRs are plotted in Figure 24 as a histogram. 
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to +40% interval about the 
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ZESK. 


37 Results For Coherent 8FSK 
The theoretical and the experimental bit error 
probability for coherently detected 8FSK  (8FSK) as a 


Bimeci Lon of the SNR inwdB is) shown Jn bigure 25. 
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REunourco 25 xsprobability of Bit Error (Theory. and 


Simulation) for Coherent 8FSK Detection. 


The relative difference between the theoretical 
probability of error for coherent 8FSK and the probability 
of error estimates obtained by the simulation is shown in 


Figure 26. 
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Figure 26. Percent difference between Theoretical and 


Simulation Results for Coherent Detection of 8FSK. 


The convergence of the estimate the bit error 
probability as a function of the number of symbols 
transmitted in the course of the simulation is shown in 
Figure 27. This allows us to estimate the minimum number of 
symbols required for a particular simulation and to verify 
that the minimum has been reached. For this particular 
Simulation, the convergence is evident at approximately 


posso symbols. 
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Figure 27. Convergence Plot for Coherent Detection of 8FSK. 


The accuracy of the simulation estimate for the bit 


error prepare Lor Coherent SFSK Is shown in Figure 2958 


The ‘oscillatory’ convergence to approximately zero as the 


number of symbols increases is evident. 
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Convergence Accuracy Verification 
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ceo ss 


The differences between the theoretical bit error 
probabilities for coherent 8FSK and the simulation estimates 
for various SNRs are plotted in Figure 29 as a histogram. 
All estimates are within the -22% to +0% interval about the 
theoretical values, with most differences between -7% and 
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4. Results For Noncoherent BFSK 
The theoretical and the experimental bit error 
probability for noncoherently detected binary FSK (BFSK or 


2FSK) as a function of the SNR in dB is shown in Figure 30. 
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Figure 30. Emobabrirt LOL Error (Theory and 
Simulation) for Noncoherent BFSK Detection. 


The relative difference between the theoretical 


probability of error for  noncoherent BFSK and the 
probability of error estimates obtained by the simulation is 


shown in Figure 31. 
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Figure 31. Percent difference between Theoretical and 


Simulation Results for Noncoherent Detection of BFSK. 


The convergence of the estimate of the bit error 
probability as a function of the number of symbols 
transmitted in the course of the simulation is shown in 
Figure 32. This allows us to estimate the minimum number of 
symbols required for a particular simulation and to verify 


that the minimum has been reached. For this particular 
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simulation, the convergence is evident at approximately 
2.3:10* symbols. 
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Figure 32. Convergence Plot (100 errors) for Noncoherent 


Detection o BESK. 


The accuracy of the simulation estimate for the bit 
error probability for noncoherent BFSK is shown in Figure 
33. In this case, the convergence is to a non-zero value 


(approximately -30$ off) as the number of symbols increases. 
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Convergence Accuracy Verification 
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Figure 33. Percent Difference Between Theory and Sim- 
ulation for Noncoherent BFSK as a function of the Num- 


ber of Symbols. 


The convergence of the estimate of the bit error 
probability as a function of the “mumber of "Symbols 
transmitted in the course of the simulation using 200 
minimum acceptable errors is shown in Figure 34. This allows 
us to estimate the minimum number of symbols required for a 
particular simulation and to verify that the minimum has 
been reached. For chars particular Simulation, the 


convergence is evident at approximately 8:10* symbols. 
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Figure 34. Convergence Plot (200 errors) for Noncoherent 


Detection of BFSK. 


The accuracy of the simulation estimate for the bit 
error probability for noncoherent BESK is shown -ing Figure 
35. In this case, the convergence is again to a non-zero 
value (approximately -30% off) as the number of symbols 


increases. 
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Convergence Accuracy Verification 
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Figure 35. Percent Difference Between Theory and Simulation 
for Noncoherent BFSK (200 Errors) as a function of the 


Number of Symbols. 


The convergence of the estimate of the bit error 
Probability as a function of thé number of symbols 
transmitted in the course of the simulation using 300 
minimum acceptable errors is shown in Figure 36. This allows 
us to estimate the minimum number of symbols required for a 


particular simulation, and to verify that the minimum has 
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been reached. For this particular simulation, the 


convergence is evident at approximately 9:10* symbols. 
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Figure 36. Convergence Plot (300 errors) for Noncoherent 


Detection of BFSK. 
The accuracy of the simulation estimate for the bit 


error probability for noncoherent BFSK is shown in Figure 


37. In this case, the convergence is again to a non-zero 
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value approximately "9G off) as the number sof symbols 


increases. 
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Figure Ey Percent Difference Between Theory and 


Simulation for Noncoherent BFSK (300 Errors) as a function 


of the Number of Symbols. 


The differences between the theoretical bit error 
probabilities for noncoherent BFSK (100 errors case) and the 
simulation estimates for various SNRs are plotted in Figure 
38 as a histogram. All estimates are within the -22% to +87% 
interval about the theoretical values, with most differences 


between -22% and 0%. 
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Figure 38. Histogram of Simulation Errors (100) for 


Noncoherent BFSK Detection. 


5. Results For Noncoherent 4FSK 
The theoretical and the experimental bit error 
probability for noncoherently detected 4FSK as a function of 


the SNR in dB is shown in Figure 39. 
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Eres 39. Probabili Nor DEO (Theory and 


Simulation) for Noncoherent 4FSK Detection. 


The relative difference between the theoretical 
probability of error for noncoherent  4FSK and TER 
probability of error estimates obtained by the simulation is 


shown in Figure 40. 
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Figure 40. Percent difference between Theoretical and 


Simulation Results for Noncoherent Detection of 4FSK. 


The convergence of the estimate of the bit error 
probability as a function of the number of symbols 
transmitted in the course of the simulation is shown in 
Figure 41. This allows us to estimate the minimum number of 
symbols required for a particular simulation and to verify 


that the minimum has been reached. For this particular 


gl 


simulation, the convergence is evident at approximately 


9:10? symbols. 
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Figure 41. Convergence Plot for Noncoherent Detection of 


4FSK. 


The accuracy of the simulation estimate for the bit 
error probability for noncoherent 4FSK is shown in Figure 
42. In this case the convergence is again to a non-zero 
value (approximately -20$ off) as the number of symbols 


increases. 
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Figure 42. Percent Difference Between Theory and 


Simulation for Noncoherent 4FSK as a function of the 


Number of Symbols. 


The differences between the theoretical bit error 
probabilities for coherent 4FSK and the simulation estimates 
for various SNRs are plotted in Figure 43 as a histogram. 
All estimates are within the 5% to 70% interval about the 
theoretical values, with most differences between 5% and 
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Figure 43. Histogram of Simulation Errors for Noncoherent 


4FSK. 


6. Results For Noncoherent 8FSK 
The theoretical and the experimental bit error 
probability for noncoherently detected 8FSK as a function of 


the SNR in dB is shown in Figure 44. 
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tion) for Noncoherent 8FSK Detection. 


The relative difference between the theoretical 
probability of error for  noncoherent  8FSK and the 
probability of error estimates obtained by the simulation is 


shown in Figure 45. 
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Figure 45. Percent difference between Theoretical and 


Simulation Results for Noncoherent Detection of 8ESK. 


The convergence of the estimate of the bit error 
probability as a function of the number of symbols 
transmitted in the course of the simulation is shown in 


Figure 46. This allows us to estimate the minimum number of 
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symbols required for a particular simulation and to verify 
that the minimum has been reached. For this particular 


simulation, the convergence is evident for approximately 


5:10? symbols. 
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Figure 46. Convergence Plot for Noncoherent Detection of 


DIU 


The accuracy of the simulation estimate for the bit error 


probability for noncoherent 8FSK is shown in Figure 47. In 


1 


this case, the convergence is again to a non-zero value 


(approximately -13% off) as the number of symbols increases. 
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Figure 47. Percent Difference Between Theory and 


Simulation for Noncoherent 8FSK as a function of the 


Number of Symbols. 


The differences between the theoretical bit error 
probabilities for coherent 4FSK and the simulation estimates 
for various SNR's are plotted in Figure 48 as a histogram. 


All estimates are within the -12$ to +5% interval about the 
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theoretical values, with most differences between -7% and 
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B. OBSERVATIONS 
Simulation results for coherent MFSK in the presence of 
AWGN agree very well with theory. As can be seen from the 


convergence plots, in all cases of coherent MFSK detection 


19 


the differences between the bit error probability estimates 
and the theoretical values converge to zero as the number of 
symbols increases, which means that the simulation estimates 
dre cobPreceEEHowevemNeturseswenot the case for noncoherent 
MFSK. The simulation results are too small for BFSK and 8FSK 
(bit error probability is underestimated) compared to the 
theoretical ones. Only in case of 4FSK are the simulation 
results higher (bit error probability is overestimated) 
than the theoretical ones. The simulation results converge 
to values that are approximately: 

e 30$ lower for noncoherent BFSK 

e 20$ higher for noncoherent 4FSK and 

* 13% lower for noncoherent 8FSK. 

The simulation results for both noncoherent 4FSK and 
8FSK are within expected bounds. Only noncoherent BFSK is 
not within expected bounds for the 95.4% confidence 
intervals for 100 errors of Monte Carlo Simulations. An 
increase in the observed number of errors to 200 and 300 
does not alter the simulation results, which still converge 
to the same somewhat incorrect values. This can be seen in 
Figure 33, 34, 35, and 36 fot. nencoherentesBESK whemer 
despite the fact that 100, 200, and 300 errors were 


observed, the bit error probability still converged to about 
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the same (incorrect) value. This suggests a systematic error 
in the Communications Toolbox for noncoherent BFSK. 

The summary of the accuracy of the simulation for both 
coherent and noncoherent MFSK detection is presented in 


Tables 4,5, and 6. 


= Bit Error Probability | Coherent BFSK | NoncoherentBFSK 
Es 4 
Max. Percent Error 9.065 87.28 
Min. Percent Error -18.209 -21.12 
LN c o S ipm 


Table 4. Comparison of the Simulation and the Theoretical 
Results for Coherent and Noncoherent BFSK Detection. 












From Table 4, we note that the standard deviation is 
much larger for the noncoherent MFSK. This means that the 
simulation results are less accurate for noncoherent BFSK 


detection. 


e Bit Error Probability | Coherent 4FSK |NonCoherent 4FSK 
Mean Percent Error 13.628 18.946 
2. | Ma 


x. Percent Error ae / 25 4 T0 







2 Min. Percent Error -10.883 69.43 


E Standard Deviation 14.109 17.931 


Table 5. Comparison of the Simulation and the Theoretical 
Results for Coherent and NonCoherent 4FSK Detection. 





SN 


From Table 5, we note that the mean values and the 
standard deviations for coherent and noncoherent 4FSK are 
closer to each other (the mean values in the absolute 
sense), so the accuracy of the simulation is comparable for 


both coherent and noncoherent 4FSK. 


- Bit Error Probability |Coherent 8FSK |Noncoherent 8FSK 


Mean Percent Error -8.207 -1.666 


Table 6. Comparison of the Simulation and the Theoretical 
Results for Coherent and Noncoherent 8FSK Detection. 










From Table 6, we note that the mean values and the 
standard deviations for coherent and noncoherent 8FSK are 
Close to each other (the mean values in the absolute sense), 
so the accuracy of the simulation is in general comparable 
for both coherent and noncoherent 8FSK. It is interesting to 
note that as the number of symbols M increases from 2 to 8 
the simulation results for noncoherent MFSK become more 
accurate while the simulation results for coherent MESK 


become less accurate. 
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V. SIMULATION AND PERFORMANCE ANALYSIS FOR CO-CHANNEL 
INTERFERENCE 
A. INTERFERENCE IN DIGITAL COMMUNICATION SYSTEMS 

Interference is a factor limiting the performance of 
digital communication systems. Co-channel interference is 
one of major types of system-generated interference and 
refers to the degradation caused by an interfering waveform 
appearing within the signal bandwidth. It can be introduced 
in a variety of ways, but most commonly by other users of 
the same portion of RF spectrum operating similar types of 
equipment [Ref 3]. In Chapter 3 MFSK digital communication 
systems were modeled, and in Chapter 4 the performance of 
these MFSK systems operating in the presence AWGN was 
verified. Since theoretical results for bit error 
probabilities are available only for the cases of AWGN, one 
must resort to simulations to estimate the bit error 
probabilities for MFSK systems with both AWGN and co-channel 
interference. 

In this chapter computer simulation results are 
provided for BFSK, 4FSK, and 8FSK communication systems 
operating in the presence of wideband, AWGN and co-channel 
interference. The results are the averages of two computer 


simulations for each case. The sampling rate was selected as 


o 


twice the minimum sampling rate (an oversampling factor of 
two), and simulations were run until at least 100 errors 
Vememonschved sin jmder@to prevent» “out of memory’ errors, 
the data sequences were limited to 10° symbols for each 
Simulation. If less than 100 errors were observed in a 
Simulation, the sequence was repeated until a sufficient 
number of errors occurred. Finally, the computer simulation 
results were obtained for a range of “AB CO tlio daB Simo 
signal-to-noise and symbol signal-to-interference (i.e., 
jamming) ratios. 
BE PERFORMANCE OF MFSK DIGITAL COMMUNICATION SYSTEMS IN 
THE PRESENCE OF AWGN AND CO-CHANNEL INTERFERENCE 

ine Order to specify the performance of MESK 
COMM c du MESvstemnscoperatrngern9bhe presence AWGN andas 
channel interference, two sets of plots are provided for M = 
2, 4, and 8, for both coherent and noncoherent detection. 
The bit error probability as a function of the signal-to- 
noise ratio (SNR) with the signal-to-interference ratio 
(SJR) as a parameter is shown in the first set, while the 
second set shows the bit error probability as a function of 
the signal-to-interference ratio (SJR) with the signal-to- 


noise ratio (SNR) as a parameter is shown in the second set. 
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qu Results For Coherent Binary FSK (BFSK) 

The probability of bit error for coherent BFSK as a 
function of the signal-to-noise ratio SNR is shown in Figure 
49. The BER increase as SNR decreases, which only happens 


for negative SJR, is discussed in Appendix E. 
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The probability of bit error for coherent BFSK as a 
function of the “aa@qnal-to-interfemence rdigies Sikes saemae in 
Figure 50. Ten curves for signal-to-noise (SNR) ratios from 
-o0 dB &e +e dB in 1.89 dB increments are shown. The solia 
line (unit wstep tunetsron)  rswthe uelbabusli ¿"OL Die Son 
for the case when only signal and interference are present 
(AWGN). Note that although the unit step has the value of 
nao positive SJR, the value of kihas keen sele red 
for SJR > 0 such that a logarithmic scale that shows the 
Simuleaemon results,well “all S920") can. bemused.” When the 
interference power is less than the signal power, (positive 
SJR) the interference has no effect, and the probability of 
error is zero because the receiver always selects the larger 
of the detected symbols. On the other hand, if the 
interference power is larger than the signal power, the 
probability of emror is 0.5 since (statistacally) halt of 
the interference symbols will be opposite to the signal 
symbols and will be selected by the receiver because of 
their larger power. Consequently, the receiver will be 
making incorrect decisions 50% of the time. 

In the noise-free case, the probability of bit errors 
caused by co-channel interference exhibits step-like 
behavior (from O to 0.5) with the threshold at SJR = 0 dB. 


As the signal-to-interference ratio increases (that is, as 
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the interference power decreases relative to the signal 
power), the probability of bit error decreases at first, but 
then becomes essentially constant (the curves become nearly 
horizontal lines). This indicates that there are roughly two 
regimes of operation: interference-dominated and noise- 
dominated. Although the curves do not change abruptly from 
decreasing monotonically to remaining constant, approximate 
values that separate the two regions may be defined by 
inspection. For example, the 7-th curve (for SNR=6.333 dB) 
becomes flat for SJR > 12 dB meaning that the performance of 
the system is interference dominated for SJR < 12 dB. 

In general, the condition SNR = SJR seems appropriate 
as the border between the interference-dominated (SNR>SJR) 
and noise-dominated (SJR>SNR), regimes of operation for 
lower values of the error probability. This is because the 
curves for the bit error probability are steep for higher 
values of SNR and SJR, meaning that a small change in SNR or 
SJR causes a large change in the bit error probability. 
Therefore, for small values of the bit error probability 
(larger values of SNR and SJR) the smaller of the two ratios 


(SNR or SJR) controls the bit error probability. 
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2. Results For Coherent 4FSK 


The probability of bit error for coherent 4FSK as a 


function of the signal-to-noise ratio SNR is shown aie euiee 


um. 
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Ten curves for symbol sSignal-to-interference (SJR) 


ratios from -5 dB to +15 dB in 2.223 dB increments are shown 
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Psa ne solixdsmne Lsethesprobabilitysof error for 
4FSK without co-channel interference (AWGN). We note the 
w matic increase in the bit error probability due tomco- 
channel interference relative to AWGN only. The probability 
of bit error for coherent 4FSK as a function of the signal- 


to-interference ratio SJR is shown in Figure 52. 
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Ten curves for symbol signadghocuolscaeeNbsNSbucsNECOn 
-5 dB &o +15 dB im 2.223 dB increments are shewn@ainm Figure 
5o The solid line (unite step function) shows the 
probability of bit error for the case when only interference 
is present (no noise). Since the plots in Figures 50 and 51 
represent the probability of bit error for coherent 4FSK, 
the ‘unit step’ may also be defined on per-bit basis. For M 


= 4 the unit step is shifted to the left by a factor of 


10-log- e which is approximately 3 dB. Also, although the 
og 


unit step has the value of zero for positive SJR, the value 
ars selected such thatla lega rita cs lema SOY 
tuac" abNtijonvcresults (which are all > 105 ) cam Dmi cd 

As was found in Section 1 (Results for Coherent BFSK), 
there are roughly two regimes of operation: interference- 
dominated and noise-dominated. Although the curves do not 
change abruptly from decreasing monotonically to being 
horizontal, approximate values that separate the two regions 
may be again identified by inspection. For example, the 8- 
th curve (for the symbol bit SNR=7.545 dB) becomes flat for 
SJR > 12 dB, meaning that the performance of the system is 
noise dominated for SJR > 12 dB. In general, as found for 
coherent BFSK, SNR = SJR is appropriate as the border 


between the interference-dominated (SNR>SJR) and noise- 
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dominated (SJR>SNR) regimes of operation for low values of 
Aa error probability. Therefore, for low values of bit 
error probability (larger values of SNR and SJR) the smaller 
of the two (SNR or SJR) by and large controls the bit error 
probability. 

S" Results For Coherent 8FSK 

The probability of bit error for coherent 8FSK as a 
function of the signal-to-noise ratio SNR is shown in Figure 
53. Ten curves for symbol signal-to-interference (SJR) 
ratios from -5 dB to +15 dB are shown. The solid line is the 
prebobility of error for the “coherent ~8FSK “without coz 
channel interference (AWGN). We note the dramatic increase 
ins. error probability due to co-channel interference 


relative to AWGN only. 
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Bit Error Probability vs SNR 
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The bit error for coherent 8FSK as a function of the 


signal-to-interference ratio SJR is shown in Figure 54. 
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Coherent 8FSK with AWGN and Interference. 


Ten curves for symbol signal-to-noise (SNR) 


ratios from 


-5 dB to +15 dB in 2.223 dB increments are shown in Figure 


54. The solid line (unit step function) is the probability 


Opine cor for the Case when only signal and interference 
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are present (noonoise). Since the plotsean Figures Se ank5 
represent the probability of bit error for coherent 8FSK, 
the ‘unit step’ may also be defined on per-bit basis. For M 


= 8 the unit step is shifted to the left by a factor of 


lon OE). which is approximately 4.71 dB. Also, although 
og 


the unit step has the value of zero for positive SJR, the 
Valiicunot (Omen se becHmeseiicceca such thare logarithmic 
scale that best shows the simulation results (which are all 
2 10 ^|. couldabewased. 

As with the two previous cases, there are roughly two 
regimes of operation:  interference-dominated and noise- 
dominated. Although the curves do not change abruptly from 
decreasing monotonically to being horizontal, approximate 
values that separate the two regions may be identified by 
inspection. For example, the 8-th curve (for SNR=5.784 dB) 
becomes flat for bit SJR > 10 dB meaning that the 
performance of the system is noise dominated for bit SJR > 
10 dB. In general, as with the two previous cases, the 
condition SNR = SJR seems to be the border between the 
interference-dominated (SNR>SJR) and noise-dominated 
(SIR>SNR) regimes of operation for lower values of bit error 
probability. Therefore, for small values of the bit error 


probability (larger values of SNR and SJR) the smaller of 


S 


MIMO SNR or SJR) by and large controls the error 
Probability. 

4. Results For Noncoherent Binary FSK (BFSK) 

mae probability of biterror for noncoherent BFSK as a 
function of the signal-to-noise ratio SNR is shown in Figure 
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Ten curves for signal-to-interference (SJR) ratios from 
-5 dB to +12 dB in 1.89 dB increments are™shown in mMotre 
oo. The solid line ds the probability omewbit “rotor 
noncoherent  BFSK without co-channel interference (AWGN 
only). We note the dramatic increase in the bit error 
probability due to co-channel interference relative to AWGN 
ORME ec 

The probability of bit error for noncoherent BFSK as a 
function of the signal-to-interference ratio SJR is shown in 
Figure 56. Ten curves for signal-to-noise (SNR) ratios from 
-5 ROB to +12 dB are shewn. -The solid line (unit step 
function) is the probability of bit error for the case when 
only interference is present (no noise). Note that the unit 
step has the value of zero for positive SJR but the value of 
10) "has been selected sof that a logarithm m- ICat 
shows the simulation results (which are all > 10” "can be 
used. As with the previous cases examined, there are roughly 
two regimes of operation: interference-dominated and noise- 
dominated. Although the curves do not change abruptly from 
decreasing monotonically to remaining constant, approximate 
values that separate the two regions may be identified by 
inspection. For example, the 7-th curve (for SNR=6.333 dB) 
becomes flat for SJR » 12 dB meaning that the performance of 


the system is noise dominated for SJR » 12 dB. In general, 


S 


psabDenMsec tae condition SNR = SJR iswethe border between the 
interference-dominated (SNR>SJR) and noise-dominated 
CJR SNE) regimes Om@Meperation for lower valweswof the error 


Probability: 


Bit Error Probability vs SJR. 
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Figure 56 Probability Ofm BIL Error versus SUR OMNE 


Noncoherent BFSK with AWGN and Interference. 
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Therefore, for small values of the error probability 
(larger values of SNR and SJR), the smaller of the two (SNR 
or SJR) by and large comtrols=the error probability. 

5. Results For Noncherent 4FSK 

The probability of bit error for noncoherent 4FSK as a 


function of the signal-to-noise ratio SNR is shown in Figure 
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haguce 57 Troba bility ofa Bit Error versus SNR. fer 


Noncoherent 4FSK with AWGN and Interference. 
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Eleven curves for symbol signal-to-interference (SJR) 
lcu cSEEcHESNdBEUOVI5SdBXn- 1-89 dB increments are shown 
in Figure 57. The solid line is the probability of bit error 
for the noncoherent 4FSK without co-channel interference 
(AWGN only). We note the dramatic increase in the bit error 
probability due to co-channel interference, relative to AWGN 
only. 

The probability of bit errřər for nomcoherent 4FSK as a 
function of the signal-to-interference ratio SJR is shown in 
Figure 58. Eleven curves for symbol signal-to-noise (SNR) 
ra ioc r om ceo +15 dB inf- 89 dB are shown. The" "Ssolrd 
linen unit step function) is the probability of bit error 
for the case when only interference is present (no noise). 
Since the plots in Figures 57 and 58 represent the 
probably” ot. b1it error for moncoherent 4FSK, the "unb 


step’ may also be defined on per-bit basis. For M = 4 the 


unit step is shifted to the left by a factor of A) 
08 


which is approximately 3 dB. Also, although the unit step 
pasen alue to zero for positive SJR, the value of 10m 
has been selected so that the logarithmic scale best shows 
MCN uU Ldcneresubts “which are all = 207). As with the 
previous cases examined, there are roughly two egies. ee 


operation: interference-dominated and noise-dominated. 


100 


0.503333, Bit Error Probability vs SJR 
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Eure" 38. Probability, of Bit Error versus roo 


Noncoherent 4FSK with AWGN and Interference. 


Although the curves do not change abruptly from 
decreasing monotonically to remaining constant, approximate 
values that separate the two regions may be identified by 


inspection. For example, the 9-th curve (for SNR-8.99 dB) 


ING 


bosm EOL AE SJR > 12. dB, meaning that the 
performance of the system is noise dominated for bit SJR > 
12 dB. In general, as with previous cases, SNR = SJR is the 
border between the interference-dominated (SNR>SJR) and 
noise-dominated (SJR>SNR) regimes of operation for lower 
values of the error probability. Therefore, for small values 
of the error probability (larger values of SNR and SJR) the 
smaller of the two (SNR or SJR) by and large controls the 


error probability: 


6. Results For Noncoherent 8FSK 

The probability of bit error for noncoherent 8FSK as a 
function of the signal-to-noise ratio SNR is shown in Figure 
59. Ten curves for signal-to-interference (SJR) ratios from 
-5 dB to +15 dB in 2.223 dB increments are shown. The solid 
line is èēhe probability of error for. noncoherent 8FSK 
without co-chanel interference (AWGN only). We note the 
dramatic increase in the bit error probability due to co- 


channel interference relative to the AWGN only. 
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Figure. 59. Probability ot MBit. EDEPOL. Versus Sb OT 


Noncoherent 8FSK with AWGN and Interference. 
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The probability of bit error for noncoherent 8FSK as a 
function of the signal-to-interference ratio SJR is shown in 


Figure 60. 


Bit Error Probability vs SJR 
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BER int only. (SJR. pitgp > €PS ) 
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Figure 60. Probability of Bit Error versus SJR fer 


Noncoherent 8FSK with AWGN and Interference. 


Ten curves for symbol signal-to-noise (SNR) ratios from 


-5 dB to +15 dB in 2.223 dB increments are shown in Figure 
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60. The solid line (step function) shows the probability of 
bit error for the cases when only interference is present 
(no noise). Since the plots in Figures 58 and 59 represent 
the probability of bit error for noncoherentWNSESEROwBhe “Unit 
step’ may also be defined on per-bit basis. For M = 8 the 


unit step is shifted to the left by a factor of 


0-Iog-( OR). which is approximately 4.71 dB. Also, although 
08 


the unit step has the value of zero for positive SJR the 
value of 10° has been selected such that a logarithmic 
scale that shows the simulation results (which are all > 
Ma usd. As  witheithe previous Casco i ex mied, 
there are roughly two regimes of operation: interference- 
dominated and noise-dominated. Although the curves do not 
change abruptly from decreasing monotonically to being 
horizontal, approximate values that separate the two regions 
may be identified by inspection. For example, the 9-th curve 
(for SNR=8.007 dB) becomes flat for bit SJR > 10 dB meaning 
that the performance of the system is noise dominated for 
bit SJR > 10 dB. In general, as before, SNR=SJR is the 
border between the interference-dominated (SNR>SJR) and 
noise-dominated (SJR>SNR) regimes of operation for lower 
values of the error probability. Therefore, for small values 


of the error probability (larger values of SNR and SJR) the 


TOS 


samademotothe two (SNR or SJR) by and large controls the 


lat error probability. 


Cc. OBSERVATIONS 

uu ccossmeurroneresults for the bak error probability as 
a function of SNR with SJR as a parameter can be used to 
determine the increase in bit error probability due to co- 
channel interference (referred here as ‘jamming’). In the 
case of intentional jamming, the objective is to increase 
the number of bit errors. In such a case, the jammer needs 
to minimize the receiver’s signal-to-interference ratio as 
much as possible. The maximum signal-to-interference ratio 
that still allows for the desired increase in the 
probability of bit errors can be determined from the 
Simulation results. Given the desired jamming-generated 
increase in bit error probability, the required jammer power 
can be calculated (also given the other relevant parameters 
such as the distance from the target receiver, the receiver 
antenna gain in the direction of the jammer, and the like 
are also known). For example, from the curves for the 
noncoherent 8FSK case, a -10 dB signal-to-interference ratio 
results in a bit error probability of 0.5, which means that 
the communications is essentially impossible. As SJR 
increases, the bint error probability decreases. 


Nevertheless, bit error probability for both noise and 


OO 


interference is always larger than the bit error probability 
for noise alone. For example, for SJR-10.23 dB and 
SNR=10.23 dB, the bit error probability is about 20 times 
larger than the bit error probability for noise alone. 

In summary, assuming that the operational SNR for the 
MFSK communication channel is known, the family of curves of 
P. versus SNR with SJR as a parameter can be used to 
determine the increase in the probability of error due to 
various levels of co-channel interference. In the case of 
intentional jamming with another synchronized MFSK signal, 
the curves may be used to determine the jammer-on-target 
requirements for a certain increase in Dit error 
probability. Subsequently, for a given scenario (jammer-to- 
target distance, etc.) the jammer requirements such as 
jammer power, antenna gain, etc. can be calculated. 

The second set of curves gives the bit error 
probability of MFSK communication systems in the presence of 
both noise and co-channel interference as a function of SJR 
and with SNR ratio as a parameter. The effect of co-channel 
interference is to introduce errors in transmission that are 
quantified by the symbol error probability Pa or the Bit 
Error Rate Pp (BER). "For "BFSK the largest bit "error 
probability is 0.5 or 50%. For MFSK the largest symbol error 


probability rs Given by: 


UNO 


Max (P.) = 1 - I/M uh 


where M is the number of symbols. The maximum bit error 


probability is given by: 


1 M 1 
M p,)2 ————(1-—)20.5 See 
ax (Pp) E "2 ( ) 


Therefore, the maximum bit error probability for MFSK is 
also 0.5, and in the noise-free case (SNR tends to 
ME) tche theoretical bit error probability due roces 
channel interference is a unit step function. The curves 
for bit error probabilities when noise and co-channel 
interference are both present tend to 0.5 for negative 
Signal-to-interference ratios (in dB) and tend 
asymptotically to zero (no errors) for large positive SJR 
ratios. The curves for negative SNR’s are essentially ‘flat’ 
Since the noise ‘controls’ the errors. On the other hand, if 
the SNR is large and positive (the noise power is small), 
the curves tend to the ‘unit step’ function of the 
interference-only case. The most realistic cases of 
comparable SNR and SJR values are those that fall in between 


the ‘flat’ curves at high (close to 0.5) values of 
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probability of error and the “unit” step Eu DNE M 


interference-only case. 








VI. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

In this thesis MFSK communication systems with 
coherent and noncoherent detection were simulated in the 
time domain using MATLAB Simulink and the Communications 
Toolbox. The bit error probabilities for coherent and 
noncoherent detection of BFSK, 4FSK, and 8FSK have been 
determined by Simulation and verified against the 
theoretical values for the case of AWGN. The MFSK models 
were modified to include the co-channel interference and the 
bit error probabilities obtained for the ranges SNR and SJR 
INaeme-oudp to TIo dB: 

Simulation results for coherent detection of  BFSK, 
4FSK, and 8FSK in the presence of AWGN show excellent 
agreement with theoretical results. On the other hand, 
Simulation results for the probability of bit error for 
noncoherent detection of BFSK, 4FSK, and 8FSK differ from 
the theoretical probabilities of bit error by approximately 
-30% for noncoherent BFSK, -20$ for noncoherent 4FSK and - 
13% for noncoherent 8FSK, indicating a ‘systematic’ error in 
the Communications Toolbox implementation of the noncoherent 
MFSK detection. Finally, for symbol signal-to-noise and 


Signal-to-interference ratios larger than 15 dB (SNRs and 


je 


Skeets dB), the required computational time can be 
excessive for 100 or so errors (days or even weeks on a 200 


MHZ PC): 


B. RECOMMENDATIONS 
The continuation of the research may include: 

e Derivation of the theoretical expression for the 
bit error probability for MFSK with AWGN and co-channel 
interference (BER as a function of SNR and SJR) for the 
following cases: 

e phase-locked co-channel interference 
(relative difference between the interference and signal 
equal 0), 

e random phase of the co-channel interference 


(relative to the signal) with uniform distribution 
between O and 2m, and 


e random phase of the co-channel interference 


with Gaussian distribution. 


e Verification of the theoretical Oe error 
probabilities for the above cases using Matlab/Simulink 


models for Monte Carlo type simulations. 


eZ 


APPENDIX A. MATLAB PROGRAM SUN PREPI.M FOR COHERENT MFSK 
CASE. 


This program prepares the data file for simulation runs 
in case of coherent detection of a MFSK communication system 


in presence of both additive Gaussian noise and co-channel 


to Wie E 

Crear 

nosccconm mM mMenu o c 
LNG Si ea 
T 


NOISE and INTEREERENCET I; 


num levels = input( Enter the number or freqtecneves sain: 
05 
Rs ca nulo ye LS) nu leve lis end 


T sym = anput ("Enter @he symbol «duration E LT]: 3955 
1f isempty(T sym), T sym = 1; end 


oversampling = input('Enter the oversampling factor [2]: '); 
i: sempesdewersampbrng), Cversampling 2; end 


MOSR = Input (Enter the MIN Sigma to Nol Sema Vegi a e: 
m 
ARES ETE ASNO NM El 


maxegub-—— wapüutüwenter tne MAXWNSMgNSNM to Nose ratio (ds 
dB]: '); 
if isempty(max SNR), max SNR = 15; end 


ihe se Omi y= =i 

minr MM OU nee im=e tie: MEN eS iGid oC mec 
enema db EE 

irf Seinen Sa, min SJR = -5; end 

noc IES pA En er cde 
mogerorlOoNsSB]: 'J); 

15 Semprey (max SIR), Max SIR — 10; end 


end 


if min SNR == max_SNR 


Nie 


mum messe = 1% 


EH 

Ains scc pue (Enter be Miúmber of values for SNR 
[10]: '); 

if isempty (num noise), num noise - 10; end 
end 


E NOS CRA 
A COSI 
num jam = 1; 
else 
pM LE on Det ot valves shor scum 
UA 
semper yam), num jam >= 10;—eng 
end 
else 
Do ley 
end 


IO SS input ('Enter EIN dE errore 
acceptabblesw00!: 7); 

crises Innes) amp cIrross pn 
cuumuEeor — Inpuly Enter the factor multiplying the 
number of errors [2]: '); 

PEG SeMpEV(CEbGt faekom error factor a7 end 
Is o lis = error a toma Saro DS; 

Max rEandint = input ("Enter the maximum size oí the random 
integer arrays [10°6): '); 


if isempty(max randint), max randint - 10^6; end 


Mile pum e Enter the file mame to save data (ne 
cs 


save sun dacal 
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APPENDIX B. MATLAB PROGRAM SUN FSK1.M FOR COHERENT MFSK CASE 
This program runs coherent MESK co-channel 

interference with additive Gaussian noise. 

clear 


load sun dati rhis losdseedlgasse. ni roundina 
delta freq = 1/T sym 


f max = (um levels-1)*1/T_ sym; 

CHI M = 0.5/(f max*oversampling) 

seeds = randinti3s, 1, 90999 

Sayeed. oec E ESA 

nozse™seed = seeds(2); 

ineert seed -NScecdse 

MET NO symbols = error fTaceOusmim errors, 
tre 


DE MIDEPCMSO 


delta SNR (mak SNR - min SNR) / (num noise -1); 
else 

AI "o 
end 


ONR - mincoNR S00: mumemnolse -—MIEISCHNENESNE 
noise var vect = T sym/(2*delta t) .* 10 2^ X-SNBRUBBONE 


Dane CcCeonly -— d 


it num jam > 1 
deltar JR max SJR min DID CE 
else 
A 
end 
SUR Ss e LO A A AE ESO ME 
interf gain vect - 10 .^ (-SJR720)5 
else 
Sue $ There is no Jamming so™the SJR 
ee Sm ee 
interf gain vect = 0; 
end 
BER = zeros(num_noise,num jam); 


total symbols = initial num symbols; 
LOL nO Ls efes MINO Es e 


noise var = nose var vect(ínoise Case); 


ES 


RETO 2c — linum Jam 


enough errors = 0; 

Inc USUS = 0; 

num symbols AOS yo > > 

rand int — Main immo symools max randint]s 


total symbols = num symbols 


WALES encudn "emtors == 1 
interf AO Ves lam. case)lo 
clear cerror nuno 
sim('fskco bm',num symbols) 
Ica cci Succ e Eo nube 
MUMS — nummcrr cH NS 
IN LL 
num symbols = iM s yoo Ss mln Ch Rotter 
seta! Lada = min([num symbols 


max randint]); 


Botal symbolsSe= total symbols gg 


num symbols; 


elselr (num err FO & Num cnn < min errors) 
num symbols SMG ee Fone 


num err) *ceil(total_ symbols/num err); 


total symbols - total symbols + 


num symbols; 


eise 
SemougBerrors le 
end 
end 
asco errors (noise case, jam Casc) m Num err 


noo MOS moss case, y amic ase > 


total) Symbols 


BER norse case jan Case) num err / 


total symbols 


end 
eval([' save ' file name '.snr SNR -ascii']); 
eval([' save ' file name '.sjr SJR -ascii']); 
eval([' save ' file name '.ber BER -ascil']); 
eval( li save ' [ile namei :-Ncr AUMber OL errors - 
ascir']); 
eval([' save ' file name '.nsy number of symbols - 
oe Dee eG 
end 
COC 
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APPENDIX C. MATLAB PROGRAM SUN PNC1I.M FOR NON-COHERENT MFSK 
CASE 


This program prepares the data file for simulation runs in 
case of non-coherent detection of a MFSK communication 
system in presence of both additive Gaussian noise and co- 


channel interference. 


oso This prepares the data file for sun runs 
clear 
select, “Pais 


Y 
OSE e 
" NOISE and INTERFERENCE '); 


nolse only — menu ( 


Düm levels —- input (Enter the mumberee s freeones om) mean. 
M 
ste DEM c O LE ES A 


T sym = input('Enter the symbol duration T [1]: *); 
IS CMA ISS IA Sy ena 


oversampling = input('Enter the oversampling factor [2]: '); 
BI seno Oe ts aman pe overscan 2. ena 


Min SNR = input( Enter the MIN Signali o Noise mnono E 
le 
LS cn (Mi un B = 9 ce 


MaS NR pue Ember thee MEX 5g dguodetogNOI Scrum 
dB]: '); 
rovcsenpbunaXxEoNPEnoxN SN Es cue 


Dc Emi i 

mins Jee input Enter tche MIN Signal to Interference 
DOO PUDE 

if isempty (min SJR) min SJR ES > 7 gene 

max ike: — AA ESA CEM o na CO ICI 
ramio[10 daB]: "y 

iS Say (axe e MARS UR = 10 Mena 
end 


if min SNR -- max SNR 


T 


DUO se ul; 
ese 

num noise = input('Enter the number of values for SNR 
[10]: *); 

1f isempty(num noise), num noise = 10; end 
enc 


ie Ol Semon eel 
if min SJR == max SJR 
num jam - 1; 
else 
num jam =sanput ("Enter the number of values for SJR 
[10]: '); 
if isempty(num jam), num jam - 10; end 


end 
else 
num jam = 1; 
end 
MINO ES Moi c emccroscuemumeun" om Or rors 
acceptable [100]: '); 
lf isempty(min errors), min errors = 100; end 


Sige) Lacon — Input mater the factor multiplying the 
ada errors M 

TEE In SEM E Or r aC COren or cameo = 2 end 

initial numwesymbolswmmerror factor*min errors; 

i RU ner the maximum Size of the random 
integer arrays [10%6]: '); 

Pisoni III dEODdgrmcemaxcrandsnt 1016; end 


A nte: che hile mame to Save dala [no 
esu, unis 


Sey Stig) orgie! 
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APPENDIX D. MATLAB PROGRAM SUN NC1.M FOR NON-COHERENT MFSK 
CASE 


This program runs coherent MFSK co-channel interference 


with additive Gaussian noise. 


clear 


load suni dnc solo cas al O EIN 


delta freq = 1/T_sym 


f max = (num levels-1)*1/T sym; 
delta t - 0.5/(f max*oversampling) 
seeds ="Handi nes, 17 LOO) 


Emonodge cog Seeds 


seeds (2); 


noise seed 
interf seed = seeds (3); 
initial num symbols = error factorem ferros, 


CIC 


IOE nüm noise 2r 


delta SNR - (max SNR - min SNR) / (num noise -1); 
else 
delta SNR = 0; 


ir 


end 


SNR — min SNR -* [0:num noise - 1]*delta SNR; 


el em cc I sym (2~delta t) * 10 .^ (-SNR/10); 


Ir nonse only -=r 
deem Jang l 
SAS A AS ESAS EA num jamne 
else 
delta SJR - 0; 
end 
SUE Min s Ju ctus ramos am EE e delta SIR, 
Mere LO at = SOR ZO) ; 
else 


SJR = - 100; $ There is no Jamming so the SJR 
jpdsNww-inrinlty 


i 
© 
~e 


interf gain vect 


end 


BER = zeros(num noise,num jam); 


total symbols = inwtial numysymbols; 


ONG ISLS Case- iii num noise 
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noise var - norse var Veco Melo So 


for jam case ~T RUNE M 


enough errors - 0; 

numser = 0; 

num_symbols = initial num symbols; 

camine = min([num symbols max randint]); 


tors vmbpoPFS num symbols 
Ware meno IM o ee 
incer rogain nter oain Weel (qebi ca eii 
ellscar crono 
sim(* tekneebm”, num symbols) 
mew srrorsterr eols "size (ercor Mier), 


num err = num err zi new errors, 


IAS e 
num symbols = num symbols*min errors, 
Once = min([num symbols 


max_randint]); 


totalTsymbols totalis ynmpo lsk 


num symbols; 
eS eee TANTS aia o acm ttt ce 


num symbols = (Minger rors > 
nüm err Gel (tocas mos ls num Err) 


Hn 


toral Symbols 


total Symbols + 
num symbols; 


else 


ll 
| 
"e 


enough Errors 
Ene 


end 


UMD SBE Peimmors (MOlsomeGacc, Jail case) hum err 


Dunber-erwswucolsinorsexcase,jam Casey 
total symbols 


BER(noise case,jam case) 
toral y meso Us 


num err / 


end 

eva MM oct le name sar toas 
eval ( [i Save tite names .S5)5 (Sue —aceri |); 
EM cM lle name '.ber BER -osemt' ]»; 


Jo IC te Mane T ner nunpermOr errors > 
Siesta sl )n- 


EUM ccu. inum DS emo e CTO LES yu csi = 
aser e 


end 


COC 
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APPENDIX E. THEORETICAL BER FOR COHERENT BFSK WITH AWGN AND 
CO-CHANNEL INTERFERENCE 

The BER as a function of SNR for negative values of SJR 
(interference power exceeds the signal power) exhibit 
somewhat counter-intuitive behaviour. Namely, the BER 
increases as the SNR increases (there are more errors for 
lower noise power levels). In order to resolve whether or 
not this is an artifact of the simulation, we have derived 
the expression for the theoretical BER for coherent BFSK 
with AWGN and co-channel interference which is assumed to be 


coherent with the signal. This result of the derivation is: 


-SJR -SJR 
] 10 20 +1 10 20 -1] 
A eec Site E 


[s 3b N 


where SJR and SNR are assumed to be in dB. 
Equation (E.1) allows us to plot BER as a surface with 
SNR and SJR as variables, as shown in Figure 60, for the 


ranges of SNR and SJR between -20dB and +20dB. 
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BER as a surface vs SNR and SJR 
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Figure 60. Probability of Bit Error versus SNR and SJR for 


Coherent BFSK with AWGN and Interference. 


The probability of bit error as a function of SNR and 
SJR shows a region of shallow local minima with the 
probability of bit error decreasing as the noise power 
increases. This occurs only for high values of interference 


and noise and only up to an SJR of O dB. For higher values 
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of SJR there is no local minimum, the decrease with SNR 
becomes monotonic. The conclusion is that noise actually 
reduces thei probability of bit error at low SIR This as 
mainly of academic interest since the reduction occurs at 
very high values of error probability (from about 0.5 to 
about 0.35). 

Equation (E.l) can also be used to generate conventional 
BER plots. These plots represent cross-sections of the 
surface illustrated in Figure 60 along planes of constant 
SNR, shown in Figure 61, or constant SJR, shown in Figure 


CZ. 


BER vs SNR for SJR as parameter 


0.55 1 





on SNR 20, 


Figure ol. Probability of Bitte error versus SNN for 


Coherent BFSK with AWGN and Interference. 
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En 5 vesctoreoUNREEPatdos from -20 dB to 0.4 dB in 4.1 

dB increments are shown in Figure 61. 
The theoretical BER for coherent BFSK as a function of 
the SJR is shown in Figure 62. Seven curves for SNR ratios 


from -20 dB to'4.5 dB in 4.1 dB increments are shown. 


Bit Error Probability vs SJR 
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Noise Variance 


Biguce 62. Probability of Bit Error versus SJR. for 


Coherent BFSK with AWGN. 
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